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Abstract. Correlations and fluctuations (the latter are directly related to the 2-particle 
correlations) is one of the important directions in analysis of heavy ion collisions. At the 
current stage of RHIC exploration, when the details matter, basically any physics question is 
addressed with help of correlation techniques. In this talk I start with a general introduction 
to the correlation and fluctuation formalism and discuss weak and strong sides of different type 
of observables. In more detail, I discuss the two-particle pt correlations/ (pt) fluctuations. In 
spite of not observing any dramatic changes in the event-by-event fluctuations with energy, 
which would indicate a possible phase transition, such correlations measurements remain an 
interesting and important subject, bringing valuable information. Lastly, I show how radial 
flow can generate characteristic azimuthal, transverse momentum and rapidity correlations, 
which could qualitatively explain many of recently observed phenomena in nuclear collisions. 



1. Introduction 

The physics of high energy heavy ion collisions attracts strong attention of the physics community 
as creation of a new type of matter, the Quark-Gluon Plasma, is expected in such collisions. 
During just a few years of the BNL RHIC operation many new phenomena has been observed, 
such as strong elliptic flow pQ, and suppression of the high transverse momentum two particle 
back-to-back correlations |2j . Note that both of the above mentioned discoveries at RHIC have 
been made with help of correlations techniques: in this case analyzing azimuthal correlations. 
Further understanding of the physics processes responsible for both phenomena, such as what 
happens with the energy of the away-side jet or whether the observed constituent quark 
scaling (31 |1J 03 in elliptic flow reflects the hadronization via quark coalescence, require further 
dedicated correlation studies. In general, as we try to understand the details of system evolution 
we will rely more and more on correlation measurements. 

Many collaborations recently have performed measurements of either multiplicity (such as net 
charge, K/tt ratio) or mean transverse momentum fluctuations H3E]. The interest was driven 
by a possibility to observe a dramatic change in fluctuation pattern with collision energy and/or 
centrality indicating the phase transition. There exist theoretical prediction both for an increase 
of fluctuations jSj as well as for decrease in fluctuations [HI E3] (the latter is due to an increase 
in number of effective degrees of freedom in the QGP). At this time many new observables were 
introduced in the analysis. Unfortunately, the various observables used by experiments almost 
render a valid comparison of the results difficult or even impossible. The reason for that is a 
persistent confusion about properties of different observables used to measure fluctuations. Due 
to the importance of inter experiment comparison of the results, in the next section I try to 



give a short summary of the question along with my own recommendations for future analyzes. 
This part is somewhat technical, and a reader more interested in physical results may skip it 
and start with section II. 



2. Correlation functions and fluctuations 

In a multi-particle process the correlations in particle production are described via two and 
many particle densities [111 I12j . In this talk I use only single and two-particle densities; those 
are normalized respectively by the mean multiplicity and the mean number of pairs in a given 
region of momentum space: 

dxp^ l \x) = (n) ; / dx\ I dx 2 p^ (x\, x 2 ) = (n(n — 1)) . (1) 
AX Jax JAX 

Particle densities in general depend on a particle three momentum, e.g. rapidity, transverse 
momentum, and azimuthal angle, some of the arguments can be integrated out; for brevity I 
use just x. From p^ l \x) and p^ 2 \x\,x 2 ) one constructs different types of correlation functions: 

C(x 1 ,x 2 ) = p {2) (x 1 ,x 2 ) - /3 (1) (xi)p (1) (x 2 ); B(x 1 ,x 2 ) = C ^ f ,X2 ) \ R(xi,x 2 )= j^~rnY7 — V 

(2) 

Roughly speaking C(x\,x 2 ) has a meaning of a distribution of correlated pairs. The correlation 
function B(x a ,Xb) corresponds to the distribution of particles b under condition that particle 
a is found at x a . Often, a particle of type a is called a "trigger" particle and particle of type 
b is called an "associated" particle. When particles a and b carry opposite charges, such as 
electric, baryon, or strangeness, B(x a ,Xb)dxb = 1, as there must be one associated particle 
somewhere in the momentum space that "balances" the charge of the trigger particle. Then 
B(x a ,Xb) is often referred to as balance function. The correlation function R(x\,x 2 ) may be 
interpreted as the "probability" that a given pair is correlated. This type of the correlation 
function is used most often, as it is easiest to measure. Defined as a ratio, it is one of the 
so-called robust quantities, which do not depend on single particle reconstruction efficiencies. 

Studying the correlations in nuclear collisions one usually is interested if those are different 
from the correlations in elementary collision. If one superimposes a few, N co u, nucleon-nucleon 
collisions together: 

p^ AA {x) = N cMP ^ NN {x) (3) 

p^ AA ( Xl ,x 2 ) = N coll p^ NN (x 1 ,x 2 ) + N coll (N coll - l)p^ NN ( Xl )p^ NN (x 2 ). (4) 

Then, the balance function, B(xi,x 2 ), does not change at all, and the signal in R(x±,x 2 ) gets 
diluted: R AA ( Xl ,x 2 ) = R NN (x 1 ,x 2 )/N coll . 

All of the above correlation functions can be calculated at fixed multiplicities. In this case 
they usually carry subscript n, e.q. C n (xi, x 2 ), and called semi-inclusive. In nuclear collision 
analyses it is not possible to fix the multiplicity of individual nucleon-nucleon collisions, and 
one should compare the results with inclusive NN correlation function, although one should 
be careful comparing very peripheral collisions selected on the basis of total multiplicity. In 
this case the approximation of correlations in individual NN collision by inclusive correlation 
function can break down. 



2.1. Fluctuations 

Many different quantities have been suggested for fluctuation analyses. Before further discussion, 
let us formulate in general terms what constitutes a good observable. It should be (i) 
sensitive to the physics under study, (ii) it should be defined at the "theoretical" level, be 
apparatus/experiment independent, (hi) have clear physical meaning, and (iv) it should not be 



limited in scope, provide new venues for further study. If most of the proposed observables do 
satisfy the first requirement, very few satisfy (ii)-(iv). It made basically impossible to compare 
many of the results, even published, one to another. Comparing different observables I argue 
that the "old fashion" correlation functions, and, similar the two particle transverse momentum 
correlations, (Ap^iAptp), are the best in many respects. 



2.1.1. Multiplicity fluctuation measures How the two, correlations and fluctuations, are related 
to each other? Let us look at the so-called reduced variance: 

, , - a l _ i , ,_v 5 .5 _ JAX dx llAX dx 2p2{xi,X 2 ) ( , 

u n = 7— r = 1 + \n) Rax, Rax - 7 1 — ? -. -. — ; — -, — r, (P) 

(n) J Ax dx 1 J Ax dx2Pi{x 1 )pi{x2) 

where the last equation shows the relation to two particle density. If pi(xi) does not vary 
much over the region AX, Rax gives an average of the correlation function R(x\,X2) over the 
momentum region used in the analysis. It follows that the reduced variance deviates from unity 
(which is interpreted as the value for statistical fluctuations) if and only if the correlation term 
R equals zero. 

Note that u depends on the product of (n) and R. Given R is a robust quantity, the mean 
event multiplicity is not, it depends on detector efficiency, track quality cuts, etc. Another feature 
of uj is that it can be different even if the underlying correlations are the same. Let us compare 
fluctuations of charged particles with fluctuations of only negative (or only positive) particles. 
Even if the correlations between particles have no any dependence on the particle charge, the 
reduced variance would show different values simply because only "half" of the particles are 
being considered. Lastly note that the reduced variance can not be independent of the size of 
the rapidity region AX as it sometimes assumed (it is possible only if R(xi,X2) oc 8{x\ — X2) 
which is not physical) as often people compare its values at different AX. 

Why one want to study ({n) R) instead of just R? The reason is that this product should stay 
constant as function of centrality in a case when AA collision is considered as a superposition of 
independent nucleon-nucleon collisions. In this case the correlation functions scales as l/N co u, 
(n) oc N co ii, and the product remains constant. The deviation from constant is easy to observe. 
Experimentally, the centrality is often obtained by measuring multiplicity, and it was wrongly 
assumed that such quantities as reduced variance do not depend on multiplicity in general. It 
would be desirable to report both quantities separately. Instead of (n) it would be better to use 
dn/dy, a quantity totally corrected for efficiency and only weakly dependent on AX. Besides, 
if R is known, one can check how it scales also other quantities, such as number of participating 
nucleons, N par t, or the number of binary collisions, or something else that might be related to 
the number of sources of particle production. 

To avoid the so called "volume fluctuations" (the effect of mixture of different impact 
parameter collisions in one event sample) the fluctuations in particle ratios are often used, e.g. 
K/ir, where K and n stand for the particle respective multiplicities in a given event. Again, 
many "details" often forgotten. First, even in the original paper [Sj, the fluctuations of particle 
ratios are "reduced" (approximated) to the following 

Wtf~W 1 W WW - W) + V) +^kk + R^-2R k .). (6) 

Note, that we started with statistically rather bad defined quantity as any ratio of two random 
numbers is (what do you do if denominator is zero?) and ended up with a much better behaved 
quantity defined via the correlation function. The question is why one would want to deal with 
ratios when even theoretical predictions are based on a really different quantity, which only 
approximately corresponds to fluctuations in ratios, namely only in the case when multiplicities 
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Figure 1. The charged particle correlation 
function #(771,772) for 771 = at various ISR 
energies Q2]. v^= 14 > 23 i and 63 GeV - 



Figure 2. (Ap^iAp^) for different 

centralities as function of the size of the 
pseudorapidity region |18j . 



of particles used in the denominator is large? It is much better to report directly the correlation 
of the type (Rkk + R-k-k — 2Rkw) [IS]- The latter is easier to measure (robust quantity) and 
can be analyzed for particles with very low mean multiplicities (for example, p and p). 



2. 2. Mean pt fluctuation measures 

Mean pt event by event fluctuation analyses are often performed by calculating the event-wise 
mean transverse momentum, an overall event average, and the variance: 
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where the first term in cr? , is defined by a single particle pt spectrum and is attributed 
to statistical fluctuations (the fluctuations in a case of independent particle production with 
the same single-particle distributions). The non-trivial fluctuations (non-statistical, sometimes 
called dynamical) are included in the second term, (Apt t iApt^) |14j . All the measures of 
fluctuations used by different experiments are related to this second term (two particle transverse 
momentum correlation (covariance)), although sometimes in a rather complicated way. In many 
cases such a relation involves tracking and analysis efficiencies, and if those are not reported 
separately, the final results become not comparable to other experiments/measurements. 

If all events in the event sample have the same multiplicity, then the equations above provide 
semi-inclusive mean transverse momentum and correspondingly (Apt iApt z). If the multiplicity 
vary within the event sample, there are two possibilities: one can average the semi-inclusive 
quantities over all events (what is used most often, I would call it event-by-event (EbyE) average) 
or calculate true inclusive quantities 
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The latter cam be written directly using two particle density 

J AX dxi J AX dx 2 p( 2 \x 1 , x 2 )Apt^Apt^ 
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Although quantitatively both definitions (EbyE and inclusive) would give almost indistinguish- 
able values (unless the multiplicity is too low and event sample include events with very different 
multiplicities), theoretically, it is easier to analyze inclusive correlations. Therefore my prefer- 
ence would be to "migrate" from the EbyE to inclusive observables. 

Note also that (Apt i Apt 2) can be easily generalized for study of correlation between particles 
of different type and for differential measurements such as Apt 1(971, cft^Aptpiift, <p2)- Similar to 
the multiplicity correlation, its centrality dependence can be tested against different hypothesis 
such scaling with rapidity density, number of participants, etc. 



3. Mean p t fluctuations/correlations at RHIC 

Recently STAR Collaboration has reported results on integrated two particle pt correlations 
for Au+Au collisions at different collisions energies ^B]- Fig. El shows the dependence of the two 
particle pt correlations on At? = 2\r] max \ the size of the pseudorapidity window taken around 
midrapidity. The results for different centralities are shown starting from the most peripheral 
(on the top, the strongest correlations) to the most central. The dependence on the size of 
the pseudorapidity region is rather weak. At higher collision energy (^/snn = 200 GeV) the 
dependence is even weaker. Under assumption that spectra dependence on transverse momentum 
and pseudorapidity factorize, according to Eq. the pseudorapidity dependence of (Ap^iApt^) 
should follow the one of correlation function R(r/i, eta-i). Then (Apt iApt 2) oc R. For a simple 
estimate one can fit the correlation function by the form i?(Ar?) oc 1 — a|Ar/|, which translates 
into R oc 1 — 4/3aY , where Y = (Ar]) max /2. The blue lines shown in Figures 1 and 2 correspond 
to the same a = 0.16 and describe the data well. 

Note that the slope of (Apt iApt 2) dependence on Y changes with centrality. It is noticeably 
smaller for central collisions. Taking into account the relation to the correlation function, it 
means that the correlations in more central collision widens in rj. The same conclusion can be 
derived directly from the differential measurements of the two particle pt correlation on (A77, 
Ac/)) presented at this conference (see Fig. 6 in ^Ej). I argue below that such correlations can 
be due to radial transverse flow. 

Fig. El shows STAR results on the centrality and incident energy dependence of 



(Apt^iAptp) / {pt) ■ Remarkably the results in this form (two particle relative transverse 
momentum correlations) exhibit almost no dependence on the collision energy. The insert 



shows the results for the most central collision where J (Apt^Aptp) / (pt) ~ 1-2% from lower 
SPS (Pb+Pb collisions, CERES Collaboration 0) to the top RHIC energy. The centrality 
dependence of (Apt^Apt^) roughly follows the expectations for a dilution of the correlations 
as more and more independent nucleon-nucleon collisions mixed up. For more detailed analysis 
of the centrality dependence note that in Eq. both numerator and denominator depend on 
centrality. Taking this into account one finds 

(Ap t! iAp t)2 ) AA ~ D Ncoll {Ap ttl Ap tj2 ) NN ; D Ncoll = — TT^n TVs ' (10) 

{N co ii - L)\n) NN + \n{n - i))nn 

The factor D takes into account the dilution of the correlations due to a mixture of particles 
from N co u uncorrelated NN collisions, and that in an individual AW collision the mean number 
of particle pairs, (n(n — 1)}nn, 011 average is larger than (n)^r^y. At ISR energies in central 
rapidity region (n(n — 1))nn ~ l-66(n)jyjy (see Fig. 0). Taking into account the dilution 
factor Dn co11 one find that the correlations increase for about 50% with centrality relative to the 
expectations based on the superposition of independent NN collisions (see Fig. El below) . 
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Figure 3. \/ (Apt^Apt^) / (pt) as a function Figure 4. J (Apt^Ap t ^} / (pt) (calculated 
of centrality for Au+Au collisions at different as non-statistical part in (p t ) event-by-event 
collision energies [TH] . fluctuations) for p+p and a + a collisions [Ej . 



4. Transverse radial flow and two-particle correlations 

At the first stage of a AA collision many individual nucleon-nucleon collision happen. Parton 
re-interactions lead to pressure build-up and the system undergoes longitudinal and transverse 
expansion. Transverse flow in the system creates strong position-momentum correlations in the 
transverse plane: further from the center axis of the system a particle is produced initially, on 
average the larger push it gets from other particles during the system evolution. As all particles 
produced in the same NN collision have initially the same spatial position in the transverse 
plane, they get on average the same push and thus become correlated. This picture leads to 
many distinctive phenomena, most of which can be studied by means of two (and many-) particle 
correlations |25j . 

The single particle spectra are affected by radial flow such that the mean transverse 
momentum is mostly sensitive to the average expansion velocity squared (pt)AA ~ (pt)NN + 
a>(v 2 }, (see Fig. El below) and to much lesser extend to the actual velocity profile (dependence 
of the expansion velocity on the radial distance from the center axis of the system) . The two- 
particle transverse momentum correlations ^3], (Ap^iAp^), measure the variance in collective 
transverse expansion velocity, and thus are more sensitive to the actual velocity profile. We 
employ a thermal model ^H] for further calculations. In this model particles are produced by 
freeze-out of the thermalized matter at temperature T, approximated by a boosted Boltzmann 
distribution. Assuming boost-invariant longitudinal expansion and freeze-out at constant proper 
time, one finds 
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where pt is the transverse flow rapidity, (fib is the boost direction, at = (pt/T) sinh(pt), and 
0t = (mt/T) cosh(pt). It also assumes a uniform matter density within a cylinder, r < R, 
and a power law transverse rapidity flow profile pt oc r n . Two particle spectrum for particles 
originating from the same NN collision in this picture can be written as 

dnp a i r t 2/n—l 



dp t ,idp tj2 



~ J dp t d<j) b pt /n 1 J(pt,i;T,pi,0 b )J(pi i2 ;r, / 9i,0 b ). (12) 
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Figure 5. (color online) Mean transverse momentum and two particle pt correlations according 
to the blast wave calculations. T = 110 MeV. 

It additionally assumes that during the expansion time (before the freeze-out) the particles 
produced originally at the same spatial position do not diffuse far one from another compared 
to the system size. The results of the numerical calculations based on the above equations are 
presented in Fig. [5] as function of (p 2 ) = (pt) 2 (4n + 4)/(2 + n) 2 . The results are shown for two 
different velocity (transverse rapidity) profiles, n = 2, and n = 0.5. One observes that indeed 
for all the particle types presented, (pt) depends very weakly on the actual profile. On opposite, 
the correlations are drastically different for two cases studied. 

In Fig. we compare our estimates with STAR preliminary data |18j on two particle pt 
correlations (taking into account the dilution factor, Eq. I10|) . We use (pt) and T parameters 
from t23J and assume N co n = N part /2 and (Apt^Apt^) / (pt) 2 = 0.011 (about 10% smaller than 
measured at ISR ^7] (Fig. 0J. It is observed that the transverse flow with n = 1 produces too 
strong correlations. 
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Figure 6. Comparison of the Blast Wave 
calculations for two different velocity profiles 
with preliminary STAR data ^B]- Relation 
(n(n — 1))nn = l-66(n) 2 VAr has been used. 



Figure 7. (color online) Two pion Acf) 
distribution as function of (p 2 ) in the Blast 
Wave model. Linear velocity profile and T = 
110 MeV have been assumed. 



While transverse flow generates in general an elongation of the pseudorapidity correlations 
(along with narrowing in azimuthal angle) it should lead to narrowing of the charge balance 
function due to the increase in mean pt as |20j : Ap z = mtsinh(Ay) ~ mtAy ~ const. 



Quantitatively, the effect is consistent with experimentally observed narrowing for about 15 - 
20% of the balance function width with centrality [21] and with centrality dependence of the net 
charge fluctuations [22]. As all particles from the same NN collisions are pushed in the same 
direction they become correlated in azimuthal space. The correlations can become really strong 
for large transverse flow as shown in Fig. (for particles originated from the same NN collision). 
Our estimates show that the azimuthal correlations generated by transverse expansion could be 
a major contributor to the non-flow azimuthal correlations |24) . 

The above described picture of AA collisions has many interesting observable effects, only 
a few mentioned here. The picture become even richer if one looks at the identified particle 
correlations. Many questions require a detailed model study, but the approach opens a 
potentially very interesting possibility to address the initial conditions and the subsequent 
evolution of the system created in an AA collision. 



5. Conclusion 

The correlation techniques are a very powerful tool in our quest to understand multiparticlc 
production in nuclear collision. They have been proven to lead to many discoveries in the past 
and promise even more to come. 
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